Excitations with fractional spin less than | in frustrated magnetoelastic chains 
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We study the magnetic excitations on top of the plateaux states recently discovered in spin-Peierls 
systems in a magnetic field. We show by means of extensive density matrix renormalization group 
(DMRG) computations and an analytic approach that one single spin-flip on top of M = 1 — ■£? 
(JV = 3,4, ...) plateau decays into N elementary excitations each carrying a fraction of the spin. 
This fractionalization goes beyond the well-known decay of one magnon into two spinons taking place 
on top of the M = plateau. Concentrating on the | plateau (N = 3) we unravel the microscopic 
structure of the domain walls which carry fractional spin-|, both from theory and numerics. These 
excitations are shown to be noninteracting and should be observable in x-ray and nuclear magnetic 
resonance experiments. 



PACS numbers: 75.10 Pq, 75.10 Jm, 75.60 Ej 

Constitutive elements of condensed systems and their 
interaction laws are all well known and in spite of this, 
modern condensed matter physics is a field where new 
fundamental concepts arise continuously, mainly due to 
strong correlation effects. The clue is that the emer- 
gent laws governing a system of many interacting bodies 
could have no direct relationship with the behavior of 
each individual member. In other words, the interaction 
processes could wash out the individual properties of the 
constituents and give rise to excitations of fundamentally 
new character.— Specifically, a new paradigm is now aris- 
ing in the field of strongly correlated electron systems, 
where the concept of Fermi liquid theory is not applicable 
any more. Collective excitations with quantum numbers 
essentially different than those of the individual electrons 
are now predicted and observed in a variety of systems. 

The earliest example arose in the 1970s in the study of 
conducting polymers as polyacetylenc. For this system, it 
was proposed that conduction was due to solitons carry- 
ing the electronic unit charge but no spin. The emergence 
of these quasiparticles carrying different quantum num- 
bers than the original constituents, was understood as 
a consequence of electron-phonon interactions.— Another 
example is provided by a two-dimensional layer of elec- 
trons in a high magnetic field. In the so-called fractional 
quantum Hall effect regime at a certain filling fraction 
corresponding to a plateau in the conductivity, the charge 
of the elementary excitations is a fraction of the electronic 
charge. The statistical properties of these quasiparticles 
are intermediate between fermionic and bosonic and they 
are termed "anyons."— 

Understanding the mechanisms through which collec- 
tive processes could produce excitations different in char- 



FIG. 1: Deconfined spinons in dimerized chain, connected 
black dots represent dimers - spin singlet combination of two 
neighboring spins, arrow is for free spin. 



acter than the original constituents of a solid state sys- 
tem, is currently under intense study. In particular, mak- 
ing specific predictions of the effects of these excitations 
on the experimental observations is a very important is- 
sue of modern condensed matter physics. 

In studying the properties of many-body systems, mag- 
netic systems have provided a fertile playground espe- 
cially for elucidating very important aspects of reduced 
dimensionality and strong correlations. When one flips 
the spin of an individual electron (say S z = — | — > S z — 
the total spin of the system changes by one unity, 
AS Z = 1, so one has created an excitation carrying spin 
S z = 1. In a three-dimensional system this excitation 
was ascribed to be carried by a bosonic particle called 
a magnon, a quantum of a spin wave. It came as a big 
surprise when Fadeev and Takhtajan identified the ele- 
mentary spin quantum number of a spin wave as S z = h 
in the one-dimensional world, calling them spinons.— One 
can "have a look" at these spinons if one introduces suffi- 
ciently strong frustration in the one-dimensional Heisen- 
berg chain so that the ground state becomes dimerized. 
Then the spinon acquires a finite gap and it is visual- 
ized free spin separating the different domains of 
dimerization as depicted in Fig. [T] 

Apart from the natural S z = | value (which could still 



2 



be ascribed to the individual electron) no other fractional 
values were observed in experiments. 

The present paper is devoted to the study of fractional 
spin excitations that go beyond the usual fractionaliza- 
tion of a magnon into two spinons discussed above. For 
example, as we discuss below, the excitations on top of 
the M = ^ plateau carry a fractional spin S z = ^ , which 
we dub "tertions." This fractionalization takes place due 
to collective effects in certain magnetoelastic systems un- 
der a strong magnetic field. These excitations should be 
observable in spin-Peierls systems like CuGeC>3 and in 
Ising antiferromagnetic chains susceptible to lattice de- 
formations. These tertions should condense as a soliton 
lattice in the ground state of a system under a magnetic 
field greater than the value corresponding to the M = | 
plateau. This would allow to directly observe these ob- 
jects in nuclear magnetic resonance (NMR), x-ray or neu- 
tron scattering experiments as has been the case with 
closing the zero magnetization gap.^ 

The lattice Hamiltonian of a frustrated spin chain cou- 
pled to frozen phonons in a magnetic field reads as^ 

i i 

+ J^Si-S^-H^S!, (1) 

i i 

H is measured in units where g[iB — 1, Si is the distortion 
of the bond between site i and K the spring constant 
and the first term corresponds to the elastic energy loss. 

J i sets the overall energy scale, and a corresponding 

distance scale. From now on, we fix J\ = K = 1 to get 
dimcnsionlcss energies and distances. 

Recently we have shown that plateaux can be present 
for magnetization values M — 1 — , with N = 2, 3, 4, ... 
being the length of the periodic cell of the ground state in 
units of the lattice constant. The actual presence of these 
plateau depends on the strength of frustration, except for 
the M = plateau (7V = 2) which is always present. The 
simplest nontrivial ones, at M = | (N = 3) and M = | 
(N = 4), have been observed clearly in numerical sim- 
ulations for moderate values of frustration Ji and spin- 
lattice coupling A\r- 

We have found that these plateaux are due to the next- 
to-leading transfer processes becoming commensurate, in 
first order of the spin-phonon interaction, and they ap- 
pear at special rational magnetization values in accor- 
dance with Ref. [||. For the M = -| plateau this corre- 
sponds to the process of transferring two particles from, 
say, the left to the right Fermi point, and for the M= j 
plateau, a process involving the transfer of three particles 
from the left to the right Fermi point. Those plateaux are 
generically less wide than the zero magnetization plateau 
which is caused by the doubling of the amplitude of the 
basic transfer process at M=0. 

Close to the zero magnetization plateau, the modula- 
tion of the lattice distortions breaks into domains corre- 
sponding to a soliton latticed Domain walls carry spin 




FIG. 2: Values of the (periodic) bosonic field for different M = 
~ ground states. Squares correspond to the three (classical) 
configurations and circles to the quantum counterpart, l 9 

S z = \ and are deconfincd. In analogy to the above 
picture our purpose is to study the excitations on top 
of the nontrivial magnetization plateaux at M = 1 — , 
iV = 3,4, ... to show that one spin- flip decays into N free 
fractional spin excitations, with spin S z — jr. To this 
end, we analyze the formation of a soliton lattice on top 
of the ^ plateau state, which is an up-up-down (uud) 
modulated structure in the frustrated antiferromagnetic 
spin chain coupled to adiabatic phonons. 

Fractionally charged excitations in the systems with 
commcnsurability 3 were studied in the early 1980s in 
onc-dimcnsional electron-phonon systems numcricallyifl 
and by bosonization.— In the case of 1/3 electron fill- 
ing these works identified elementary excitations car- 
rying charge and spin values (in addition to polaronic 
excitation with ordinary electronic quantum numbers): 
AQ = ±e/3, AS* = 1/2 and AQ = ±2/3, AS = 0, re- 
spectively. Our case corresponds to 1/2 electron filling, 
with completely frozen charge fluctuations. As we will 
show in this case spin excitations will be fractionalizcd 
in the units of 1/3. Since the charge field is suppressed 
there is no direct analogy between the quantum numbers 
of the excitations for electronic systems and our mag- 
netic system which is equivalent to the system of spinlcss 
electrons. 

We developed a self-consistent harmonic approxima- 
tion (SCHA) in analogy to the zero magnetization case 
and our findings are fully confirmed by extensive DMRG 
computations. 

The bosonized version of Eq. JT]) reads like 

H = Hq + H p h + H sp , (2) 

where Hq is a Gaussian part, H p h is the adiabatic phonon 
part, and H sp is the spin-phonon interaction term which, 
around M = ^ is given by£ 

- A x J dxS(x)[(3 : cos(V2?r0) : +7 : cos(2\/27r0) :], (3) 

where S(x) is the smooth part of the displacement field 
in the continuum limit and columns : • ■ • : indicate nor- 
mal ordering of the vertex operators with respect to 
the ground state with magnetization M. For kp = ^ 
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there are three incquivalent minima which are degener- 
ate, which correspond to the three different uud arrange- 
ments. In terms of the phonon and bosonic fields, they 
correspond to S p (x) = So cos[2fc^(a; + p)] (p = 0,1,2) and 
V2n(j) = 0,±^-, respectively. These three structures are 
clearly observed in the numerical simulations. 

An interesting observation is that singlets can always 
appear in domain walls because when tunnelling from the 
first vacuum to the second or third, the field rests on the 
intermediate pscudominimum in between (which turns 
out to be a portion of the quantum plateau indicated by 
yellow circles in Fig. [3]). 

Let us now analyze the excitations on top of the M = | 
plateau. In the pure spin case, the potential energy is 
given simply by V[4>\ oc J dx cos (3-v/27r0) , which also 
has three degenerate minimal. From this potential, one 
immediately concludes that the excitations on top of the 
plateau correspond to massive kinks (whenever V[</>] is 
a relevant perturbation) interpolating between these in- 
equivalent minima, and carry fractionalized spin-iJ^ 

In the spin-phonon case, the situation is more subtle, 
since now the three minima correspond to combined mag- 
netoelastic configurations as we discussed above. To see 
how fractional spin kinks arise in that case, we resort to 
a SCHA along the lines of Refs. 0311] 

Following Refs. [l3lfl4lfl5l we split <j> into classical and 
quantum components, <p — <p c + <f> q . Using the value of 
kp for M =| and keeping only commensurate terms, we 
arrive at the following potential for the classical bosonic 
field: 



dx cos 



(3\/2^ c ) 



(4) 



which led us to conclude that kinks arc similar to those 
in the pure spin case, though now both spin and phonon 
modulations must combine appropriately. Below we find 
the explicit expression for the local magnetization and 
bond modulations and compare them with our numerical 
results obtained by DMRG. 

Let us start discussing how (j> c evolves as we walk 
around a chain with periodic boundary conditions 
(PBC). We start from the vacuum corresponding to 
y/2n(j> c = 0, then we have a tunnelling of \f2Tt(j) c from 
— > =^-, at the position of the first domain wall (let 
us call this point x±), then a tunnelling process from 
4^ — > ^ takes place (at X2) and at the position of the 
third domain wall (x^) the initial vacuum is restored by 
tunnelling 4^— >0. An analytic expression for <j> c can be 
built up as a product of three soliton solutions of the 
sine-Gordon modet^ centered at xi, xi and X3. It reads 



3V2n(p c (x) = —2 



4arctan{exp[(a; — xi)/£]} 
x (2-7T + 4 arctan{exp[(x — £2 )/£]}) 
x(47r + 4arctan{exp[(x — x 3 )/^]}) (5) 

with £ being the soliton width. From Eq.© and the 
bosonization formulas connecting <p(x = ia) with S? 



fRefllTT) we extract the local magnetization of every three 
sites, 



<5*(3a?)> = -^d x c (x)-B lC os(V2^Mx) + Y a 



Sir \ 1 

S 2 cos ( 2V2^(j) c (x) + —a] + -. (6) 



As anticipated, singlets indeed appear within domain 
walls. This is because when tunnelling from one vac- 
uum to another, the field passes through the intermediate 
pscudominimum in between, which is exactly a portion 
of the quantum platea u 18 ' 19 (sec Fig. [2]). Here we would 
like to note that in the absence of the spin phonon cou- 
pling the 1/3 magnetization plaetau in the J\ — J 2 model 
was for the first time identified for stronger values of J 2 
by Okunishi and Tonegawa^S who also identified simi- 
lar fractionalized spin- 1/3 excitations around it^i They 
connected this excitation with a domain wall in the Ising 
limit, microscopically different from a singlet-core excita- 
tion that is realized in our case of additional spin phonon 
interaction and SU(2) symmetric spin exchange. 

We now undertake a numerical analysis of the lat- 
tice deformations and the local magnetization ((S*)) 
around the plateau at M = |. We have used an iterative 
method based on a DMRG procedure to solve the adi- 
abatic equation corresponding to Hamiltonian ([T]) along 
the lines stated in Ref|7|. To compare with the previous 
analytical study, we consider PBC, and the calculations 
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FIG. 3: DMRG solution for N s = 138, J 2 = 0.5, and Ai = 0.6. 
We represent with circles (squares) symbols Si (< Sf >), using 
different colors for each of the three sublattices. (a) Results 
for M = I in our system, (b) and (c) show the results for Solid 
lines in (b) and (c) correspond to the modulations obtained 
within the SCHA (RefQ). Three well-defined excitations are 
seen as what we call tertions. Open symbols and dotted lines 
correspond to DMRG and bosonization results for a second 
pattern where the central tertion is shifted. Both patterns 
have the same energy, showing that tertions are noninteract- 
ing. 
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FIG. 4: DMRG results for the same set of parameter used in 
Fig. [3] for the local spin-spin correlation function. 

With M = J 9% , we want to study the states for 
M = | and one unit of magnetization above it. The it- 
erative procedure for the Si, takes around 100 iterations 
to achieve convergence in a particular N s . Note however 
that a periodic pattern with a wavelength A = n (n 
integer) is expected for each magnetization M £ There- 
fore, to reduce the CPU time we impose such a peri- 
odicity for both M, on the Si pattern in our numerical 
calculation, choosing N s = 138. Then we study for val- 
ues of S+fot = 23 and 24, and this enforcement help us 
to obtain very accurate results for the states we are in- 
terested in. In Fig. [3] we show the results of Si and S z 
for a particular set of parameters where the plateau at 
M = \ is present. In Fig. Ifa), for Sf ot = 23, A = 3 the 
up-up-down structure is clearly seen, corresponding to a 
weak-weak-strong structure for the bonds. Figures [3](b) 
and|3fc) show magnetization and distortion respectively, 
forSf ot = 24,A = 46. 

The patterns obtained for Si and Sf are oscillatory on 
the scale of the lattice constant. We separate the lat- 
tice in three different sublattices to extract the smooth 
variations of the relevant quantities. Three different ex- 
citations are clearly identified which are characterized as 
domain walls of the uud order. As the total spin of this 
state is S z = 1 above the M = i state, each excitation 
carries S z = | and for this reason we term them ter- 
tions. Moreover, a very accurate fitting could be found 
between the DMRG results and the analytic form given 
in Eq. ([6]). Lines on Fig. [3^c) were obtained from this 
expression with parameters B\ = 0.35, B2 = 0.03, and the 
soliton width in units of three lattice sites, £ = 4.5. 

In Figs. EDJb) and[«2c), we added further DMRG results 
and the analytical fitting, now shifting the position of the 
second domain wall, and running the code again without 
forcing the periodicity. The overall coincidence between 
the bozonization and DMRG results and the fact that 
both states have the same energy confirm that the exci- 
tations correspond to noninteracting solitons with frac- 



tionalizcd spin S z = §. DMRG results for other lattice 
sizes not shown here, lend further support to this conclu- 
sion of independence, in particular since Eq. ^ perfectly 
fits in all cases the numerical results using the same set 
of constants B\,B2, and £. We also checked that exci- 
tations behave similarly for different sets of parameters 
where the M=^ plateau is present. 

Finally, let us analyze the internal structure of these 
tertions. Looking at the tertion placed at the center of 
the lattice, it can be seen in Fig. [3£c) that S z has greater 
value at position 63, and almost vanishes at sites 61-62 
and 64-65. This fact points towards singlet formation as 
we have predicted theoretically. In fact, we have calcu- 
lated the spin-spin correlation (S.; • Sj+i) shown in Fig. [H 
and we obtain that the value ~— | at the bonds around 
each tertion is centered. Depending on the system size, 
the quantum plateau portion can be longer or shorter. 

In conclusion we have shown that plateaux in magne- 
toelastic systems, independently of the mechanism that 
produce them, involve the development of a soliton lat- 
tice at the threshold. We have also shown that solitons- 
domain walls carry fractional spin values which are gener- 
ically smaller than |, in particular for the excitations 
around the M = ^ plateau, non-interacting quasiparticles 
with fractional spin S z — i arise. We have also identi- 
fied the core of the domain wall as singlets in the case 
of the M = i plateau (Fig. [5|). We hope that our pre- 
dictions will stimulate further high field experiments on 
spin-Peierls compounds. Like for the case near M = 
in the spin-Peierls material CuGe03j£ the lattice defor- 
mation we predict could be measured in x-ray or neutron 
scattering experiments. The local magnetic texture could 
otherwise be seen in NMR experiments as in Ref. @. 
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FIG. 5: Scheme of the magnetic structure of the state at 
M — i (upper panel) and a tertion (lower panel) as given by 
DMRG calculation for the same parameters of Fig. [3] The 
length of the arrows is proportional to < Sf >. 
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